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ABSTRACT 

First evidences of IR excess and disk mass accretion (strong Ha emission) around brown dwarfs 
seem to indicate the existence of circumstellar disks around these sub-stellar objects. Nothing 
is known at the present time about outflows which potentially might be launched from brown 
dwarfs, although jets are typically associated with the accretion in standard T Tauri star disks. In 
this paper we calculate the Ha emission of internal working surfaces produced by a radiative jet 
in a neutral and in a photoionizcd environment as a function of the jet parameters (the ejection 
velocity Vj , shock velocity v s , mass loss rate M and radius rj of the jet) and we provide estimates 
of the Ha luminosity for the parameters of "standard" Herbig-Haro (HH) jets from T Tauri stars 
and for the parameters expected for jets from BDs. Interestingly, we find that while the mass 
loss rates associated with jets from BDs are found to be two orders of magnitude lower than the 
mass loss rates associated with "standard" HH jets (from T Tauri stars), their velocities are likely 
to be similar. Based on our calculations, we discuss the conditions in which jets from BDs can 
be detected and we conclude that the Ha: luminosities of internal working surfaces of jets from 
BDs in a photoionized environment should have only one order of magnitude lower than the Ha 
luminosities of T Tauri jets in a neutral environment. 

Subject headings: ISM: jets and outflows - ISM: Herbig-Haro objects - stars: low-mass, brown dwarfs 



1. Introduction 

Following the discovery (Nakajima et al. 1995) 
of the first brown dwarf (BD) , several others have 
been detected in the past few years (Martin et 
al. 1997; Kirkpatrick et al. 1999, 2001; Gizis 
et al. 2001; Burgasser et al. 2003; Gizis et al. 
2003). Many of the properties of BDs are still un- 
known, and even the definition of a brown dwarf 
remains a subject of discussion. One of the fun- 
damental open questions related to BDs is their 
origin. Are they the result of a molecular cloud 
collapse followed by the accretion of circumstellar 
disks formed around a central nucleus (Elmegreen 
1999) or are they stellar embryos that, due to dy- 
namical interactions, are ejected from a multiple 
system before having time to accrete sufficient ma- 



terial from a circumstellar disk to reach the hydro- 
gen burning limit (Reipurth & Clarke 2001; Bate, 
Bonnel & Bromm 2003) ? 

It is even more difficult to define which are the 
observational evidences that could be used to dis- 
cern between the two models. At the present time, 
there are no elements that can exclude one of these 
models, and it is quite important to collect ob- 
servational evidence and/or theoretical arguments 
that can support or reject one or more of these 
models. It is of course also possible that more 
than one of the proposed mechanisms might be 
active in BD formation. 

What about the existence of jets from brown 
dwarfs ? The existence of disks is a necessary con- 
dition for the existence of jets. Evidence of accre- 
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tion disks around BDs was recently obtained by IR 
excess estimations (IR excesses being the typical 
signature of the presence of disks around T Tauri 
stars). Near-IR (3.8 pm) detections in IC348 and 
Taurus (Liu ct al. 2003), mid-IR (7 - 14 /xm) de- 
tections in p Oph and Chamaeleon (Comeron et 
al. 1998; Comeron ct al. 2000; Natta & Testi 
2001; Testi et al. 2001; Natta et al. 2002) and 
sub-millimeter and millimeter detections (Andre 
& Montmcrle 1994; Carpenter 2002; Klein et al. 
2003) started to identify a rich sample of BDs 
with circumstellar disks. Besides this, signatures 
of mass outflows are also obtained by the detec- 
tion of strong Ha emission principally in young 
BDs (with ages of a few Myr) in the a Ori region 
(Barrado y Navascues et al. 2002; Zapatero Oso- 
rio et al. 2002), Taurus (Martin ct al. 2001), and 
the TW Hydrae association (Gizis 2002). 

It is more difficult to find a straight correlation 
between the presence of jets and the formation 
mechanism of BDs. In the scenario in which BDs 
form from a direct collapse of a molecular cloud, 
they should have an accretion disk similar to the 
one around a T Tauri star, and therefore the pro- 
duction of jets would be expected. In the scenario 
in which BDs are stellar embryos ejected from a 
young, multiple system, the BDs might still have 
low mass disks around them, which potentially 
could also eject bipolar jet systems. Due to the 
fact that the possible outflow ejection mechanisms 
associated with these different BD formation mod- 
els have not yet been studied in any detail, it is not 
possible to say what a discovery of jets from BDs 
would imply about the formation mechanism(s) of 
BDs. 

A paper that is particularly relevant in the 
context of jets from BDs is the one of Quillcn 
& Trilling (1998), who showed that even proto- 
Jupitcrs might be expected to produce jets. There 
is also the work of Wolk & Beck (1990), who ex- 
trapolated the outflow luminosity vs. source mass 
correlation observed for jets from T Tauri (and 
more massive) stars to low stellar masses, to ob- 
tain predictions of the luminosities expected for 
jets from BDs. 

Which is the typical velocity of ejection of BD 
outflows ? Is the hydrodynamic process producing 
shock emission in BD jets similar to the one of 
HH jets from T Tauri stars ? Which is the typical 
luminosity of their emission ? Can we hope to 



detect BD jets ? What would be the best way to 
search for them ? 

In the present paper, we discuss a simple, the- 
oretical "internal working surface" model, from 
which we obtain predictions of the Ha luminosity 
of an HH jet travelling within a neutral medium 
or within a photoionized region (section 2). We 
use this model to obtain predictions of the Ha lu- 
minosity of jets from BDs (section 3). Finally, we 
use our results to propose observational strategies 
for detecting jets ejected by BDs (section 4). 

2. An internal working surface 
2.1. The basic geometry 

We model the knots along HH jets as "internal 
working surfaces" resulting from a time- variability 
of the ejection velocity (Raga et al. (1990), Raga 
& Kofman (1992), Kofman & Raga (1992)). 

In order to obtain analytic estimates of the 
emission expected from an internal working sur- 
face (IWS), we consider the most simple possible 
analytic model. The schematic diagram of Fig.l 
shows the flow structure of an IWS when seen from 
a frame of reference moving with the velocity v ws 
of the working surface. In this frame of reference, 
the material from the upstream and downstream 
continuous jet segments converges onto the work- 
ing surface, producing a two-shock structure. 

The most simple situation is found for the case 
in which the upstream and the downstream jet seg- 
ments have the same density pj. In this case, the 
material in both beam segments converges onto 
the working surface with the same velocity v s (see, 
e. g., Raga & Kofman 1992). 

If the two shocks are strong and radiative, the 
shock jump conditions can be written as : 

2 2 / 1 \ 

PjV s = p s C s , (1) 

where pj is the density of the continuous beam seg- 
ments (see Fig.l), v s is the shock velocity, and p s 
and c s are, respectively, the density and isother- 
mal sound speed (« 10 km s _1 ) at the end of the 
postshock cooling regions. 

Following Falle & Raga (1993), we can now esti- 
mate the separation d between the two IWS shocks 
(see Fig.l) by assuming that the material escapes 
sideways from the jet beam with a velocity equal to 
c s , and then write the inflowing/outflowing mass 
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balance of the IWS 



5-6) we finally obtain the Ha luminosity 



M in = Mom , 



(2) 



where M in = 2nrj 2 pjV s is the rate of mass com- 
ing into the IWS through the two shocks, and 
M OMt = 2irrjdp s c s , where rj is the jet radius (see 
Fig.l). From equation (2) we then obtain that the 
separation d between the two IWS shocks is given 
by: 

(3) 



2.2. Internal working surfaces in a neutral 
environment 

Let us now estimate the Ha luminosity of an 
IWS of a jet launched in a neutral region. For 
computing a simple estimate of the emission, we 
follow Raga & Kofman (1992), and assume that 
the emission comes mostly from the two shocks 
within the jet beam, and that the emission from 
the bow shock (which is pushed into the surround- 
ing environment, see Fig.l) is negligible. 

The total Ha luminosity of the IWS is given 
by: 

L [ h1 = ^ rj 2 a Ha (4) 

where an a is the luminosity per unit area (of one 
of the shocks), rj is the jet radius, and the contri- 
bution from the two IWS shocks has been consid- 
ered. 

If one considers the Ha emission from the 
"self consistent preionization" plane-parallel shock 
models of Hartigan, Hartmann & Raymond 
(1987), one finds that for the shock velocity range 
v s = 30 — > 100 km s _1 the Ha luminosity per 
unit area c jj a of the shocks can be fitted with the 
interpolation formula 

a^^xlO^ergcms- 1 ( woS^ Y 

(5) 

where v s is the shock velocity and rij (« pj/1.3niH 
in our case, see Fig.l) is the preshock number den- 
sity. 

We now write the average mass loss rate from 
the time-dependent source as 

M = irrfpjVj , (6) 

where Vj is an appropriately defined average value 
of the ejection velocity. Combining equations (4- 



l h1 = 3-0 x 10- 3 L W 
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100 kms" 
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normalized to typical parameters for a "normal" 
HH flow from a T Tauri star. 

2.3. Internal working surfaces in a pho- 
toionized environment 

Let us now estimate the Ha luminosity of an 
IWS of a jet immersed in a photoionized region. 
As has been discussed by Raga et al. (2000) and 
by Lopez-Martin et al. (2001), in many cases HH 
jets can be completely photoionized by the ioniz- 
ing radiation field from the H II region source. We 
will therefore assume that we have an IWS which 
is fully photoionized by the external, ionizing field. 

We will then compute an estimate of the Ha 
emission coming from the volume of gas in between 
the two IWS beam shocks (see Fig.l) as : 



L (P) - U 2 



"Ha 



(nrjd)n s a Ha {T a )hv Ho 



(8) 



where the first factor in parentheses on the right 
is the volume of emitting material (see Fig.l), 
n s = p s /(1.3mj) is the post-cooling region num- 
ber density, aH a (T s ) is the effective Ha recombi- 
nation coefficient (calculated at the T s « 10 4 K 
temperature of the post-cooling region photoion- 
ized material) and hvHu is the energy of the Ha 
transition. Using equations (1), (3) and (6) and 
setting c s w 10 km s _1 we then obtain 



4 P 2 = 1.68L 



M 



lO- 7 M yr- 



100 kms 



( V -I V 

VlOOkms" 1 / 



100 AU 



.lOOkms- 1 1 ; : ' 

normalized to typical parameters for a "normal" 
HH flow from a T Tauri star. 

Actually, to this Ha luminosity we should add 
the contribution from the cooling regions right be- 
hind the two IWS shocks. However, from a com- 
parison of equations (9) and (7) we can see that 
this contribution has to be negligible, at least for 
the parameters of a "normal" HH jet. 
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Another interesting feature of fully photoion- 
ized jets is that the emission from the jet beam 
itself can be quite substantial. It is straightfor- 
ward to show that for a constant density jet the 
Ha luminosity from the whole of the jet beam is 
of the order of L H ^ (the luminosity of the work- 
ing surface, see equation 9) multiplied by a fac- 
tor (L/rj)(c s /v s ) 3 (where L is the total length of 
the jet beam and c s « 10 km s _1 is the sound 
speed of the photoionized gas). For example, for 
L/rj ~ 100 and v s ~ 100 km s _1 , we would have 
a jet beam Ha luminosity of 10 % of the working 
surface luminosity. Though being quite diffuse, 
this emission from the jet beam might also be de- 
tectable. 

3. Optical emission of outflows from brown 
dwarfs 

In this section we use the results obtained in 
sections 2.2 and 2.3 to try to compare the emission 
properties of typical HH jets with the properties 
of jets ejected by brown dwarfs. 

All of the jet (or wind) ejection mechanisms 
have in common the fact that the ejection velocity 
Vj is of the order of the escape velocity v es (vj ~ 
v es ) from the surface of the outflow source. For an 
outflow from a T Tauri star we then have that: 



V es ,TT 



/2GM TT 



(10) 



where Mtt is the mass and ttt the radius of a 
T Tauri star and G is the gravitational constant. 
For a BD, we will have an outflow with a velocity 
of the order of the escape velocity 



V es ,BD — V eSy TT ■ 



M B p ■ ttt \ 
M T t ■ r B D ) 



0.5 



(11) 



where Mbd is the mass and tbd the radius of a 
BD. 

Which are the typical values for Mbd and tbd 
? Under the assumption that the BDs which 
are most likely to produce observable jets are the 
youngest and most massive ones we can reason- 
ably estimate a range of M B d ~ 0.03 - 0.075 M Q 
for the mass. It is more difficult to identify a typ- 
ical value for tbd- From cool, sub-stellar object 
atmosphere models (Baraffe et al. 1998; Baraffc 
et al. 2003; Burrows et al. 1997; Burrows et al. 



2001) one sees that for old objects (older than ~ 
1 Gyr) the radius of brown dwarfs tbd converges 
to a typical value of ~ 0.1 Rq. For younger BDs 
(of a few Myr), tbd varies within the 0.1 — 0.4 Rq 
range. 

If we then consider typical values ttt ~ 4 Rq, 
M T t = 1M for a T Tauri star, and M B d = 
0.06M Q , tbd ~ 0.3 Rq for a BD, from equation 
(11) wc obtain that u es ,BD ~ «es,TT- Therefore, 
using typical masses and radii, we find that the es- 
cape velocities from the surfaces of BDs have val- 
ues that are similar to the ones of T Tauri stars. 
We conclude that the velocities vj^tt and vj.bd 
of ejection from T Tauri stars and BDs are also 
comparable. 

Another important parameter is of course the 
shock velocity v s of the shocks associated with the 
IWS (see, equations 7 and 9). This shock velocity 
is not only important for determining the Ha lu- 
minosity of the knots, but it also determines the 
excitation/ionization of the emitted spectrum. In 
models of outflows from time-dependent sources 
the shock velocity v s has a maximum possible 
value v s ~ Awj/2 (where Avj is the full ampli- 
tude of the ejection velocity variability). 

As no clear model for the production of ejection 
variabilities has yet been studied in sufficient de- 
tail, it is not possible to find a theoretical relation 
between the ejection velocity Vj and the ampli- 
tude Avj of the ejection velocity time-variability. 
However, from observations of HH jets (see, e. g., 
the compilation of Raga, Bohm & Canto 1996) 
it is clear that the shock velocities of the shocks 
associated with these objects range from v s ^ m i n 
~ 20 km s _1 (this lower limit probably being due 
to the fact that lower velocity shocks produce very 
faint line emission and are probably not detected) 
up to Vs^max ~ 100 km s . Therefore, the am- 
plitude of the velocity variability ranges from low 
values up to Av jimax ~ 2v Symax ~ 200 km s _1 . 
In other words, the amplitude of the velocity vari- 
ability reaches values similar to the mean velocity 
vj of the outflow. It appears thus to be reasonable 
to assume that this result might also be applica- 
ble for jets ejected by BDs rather than by T Tauri 
stars. 

Therefore, jets from brown dwarfs have prob- 
ably knots with emission spectra of similar exci- 
tation/ionization to the ones of jets from T Tauri 
jets. Also, the Ha luminosities of the knots will 
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be given by equations (7) and (9) with values of Vj 
and v s similar to the ones given for T Tauri jets. 

However, the Ha luminosity is strongly affected 
by the fact that the accretion rates onto young 
brown dwarfs appears to be roughly two orders 
of magnitude lower than the accretion rates onto 
T Tauri stars (Muzerolle et al. 2003; Barrado y 
Navascues et al. 2004). One can roughly consider 
a range of M TT ~ 10~ 8 - 10~ 7 M Q yr" 1 for T- 
Tauri stars and a range of Mbd ~ 10~ 10 — 10~ 9 
M Q yr" 1 for brown dwarfs. As the mass loss rate 
in an outflow appears to be proportional to the 
mass accretion rate (Cabrit et al. 1990), this result 
directly implies that the mass loss rate of a brown 
dwarf jet is at least two orders of magnitude lower 
than the mass loss rate of a typical T Tauri jet. 

This lower mass loss rate implies that the Ha 
luminosity of a jet from a BD embedded in a 
neutral region will be two orders of magnitude 
lower than the Ha luminosity of a typical jet 
from a T Tauri star (see equation 7). Taking 
values of Vj~100 km s _1 , v s =100 km s _1 and 
M=10 -9 M yr _1 the knots along a brown dwarf 
jet in a neutral region should have Ha luminosities 
L ( ^(BD) - 3 x 10~ 5 L Q (see equation 7). 

It is not possible to make such a firm state- 
ment for the case of a BD jet embedded in a pho- 
toionized region. The first thing that we see from 
equation (9) is that a drop of two orders of magni- 
tude in the mass loss rate results in a drop of four 
orders of magnitude in the Ha luminosity (when 
going from a T Tauri to a BD jet). If we assume 
that the radius of a BD jet has a value similar to 
the one of a T Tauri jet, we would conclude that 
the knots along a photoionized BD jet should have 
an Ha luminosity L { p a (BD) ~ 1.7 x 10~ 4 L Q (see 
equation (9)). 

4. Conclusions 

We have obtained simple, analytic estimates of 
the Ha luminosities of internal working surfaces of 
jets in a neutral environment (see section 2.2) and 
of fully photoionized jets (see section 2.3). 

We have argued that the dynamical character- 
istics (i. e., jet velocity and amplitude of the ejec- 
tion velocity variability) are likely to be similar 
for jets from T Tauri stars and jets from BDs (see 
section 3). Therefore, the lower observed mass ac- 
cretion rates (which in its turn determine the mass 



outflow rates) directly imply that "neutral" brown 
dwarf jets should have Ha luminosities two orders 
of magnitude lower than the luminosities of "nor- 
mal" HH jets from T Tauri stars. This result is 
in agreement with the Ha luminosities calculated 
for BD jets by Wolk & Beck (1990), who extended 
to low source masses the Lu a vs. M source correla- 
tion found for outflows ejected from T Tauri (and 
more massive) stars. 

We find that the observed correlation can be 
interpreted in terms of our internal working sur- 
face model as a sequence of outflows with the same 
jet velocity Vj and variability amplitude Avj, but 
with different mass loss rates Mj. 

We find that the knots along photoionized BD 
jets have Ha luminosities L^(BD) ~ 10~ 4 L Q 
(see section 4). Even though such luminosities are 
4 orders of magnitude lower than the ones of pho- 
toionized T Tauri jets (see equation 9), they are 
only lower by an order of magnitude with respect 
to the luminosities of "neutral" T Tauri jets (see 
equation 7). 

This result indicates that favourable conditions 
to carry out searches for BD jets can be found 
in regions with young, low mass stars embed- 
ded within an H II region, such as M 42 (Bally 
& Reipurth 2001), the Pelican Nebula (Bally 
& Reipurth 2003) or the IC 1396N H II region 
(Reipurth et al. 2003). However, as is clear from 
the published images, in such regions it is difficult 
to detect jets against the bright and inhomogc- 
neous background of the H II region (which many 
times has ridges or bright edges which might or 
might not be jets). This problem could be over- 
come by searching for photoionized BD jets with 
a large field Fabry-Perot interferometer, in which 
even faint jets should clearly stand out from the 
H II region emission at high radial velocities. 

As a final comment, let us point out that in 
order to fully photoionizc an HH jet, it is neces- 
sary to have a relatively strong and nearby source 
of photoionizing radiation. For producing an Ha 
luminosity of ~ 1 Lq (as predicted for a jet from 
a T Tauri star, see equation 9), it is necessary 
to have ~ 5 x 10 45 s" 1 ionizing photons. In or- 
der for a working surface to intercept this number 
of photons, it is necessary to have something like 
an 05 star (producing S» « 5 x 10 49 s _1 ionizing 
photons) within a distance of ~ 50 rj or an 09 
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star (with 5» w 2 x 10 48 s _1 ) within a distance of 
~ 20 rj. The same stars could fully photoionize 
the much more puny BD jets at distances a factor 
of ~ 100 times larger than the ones quoted above. 
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Fig. 1. — Internal working surface (IWS) result- 
ing from a variability of the ejection velocity. The 
diagram shows the flow as seen from a reference 
system moving with the IWS. In this reference 
system, the material in the jet beam converges 
onto the IWS. If the up- and down-stream jet 
beam segments have the same density pj, the 
two IWS shocks have the same shock velocity v s . 
The shocked material within the working surface 
(which has a density p s and sound speed c s ) es- 
capes sonically through the sides and interacts 
with the surrounding environment, forming a bow 
shock. The radius of the jet is rj , and the separa- 
tion between the two IWS shocks is d. 
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